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initialization  errors  and  inadequate  resolution  caused  pressure  peaks  to  appear  too  low  by  about  20%. 

In  the  Mach  reflection  region  the  peak  pressures  are  in  good  agreement  with  experimental  data. 
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TRANSITION  TO  DOUBLE  MACH  STEM  FOR  NUCLEAR  EXPLOSION  AT 
104  FT  HEIGHT  OF  BURST 


I.  INTRODUCTION 

In  a  nuclear  height-of-burst  (HOB)  detonation  the  spherical  blast 
wave  reflects  from  the  ground,  initially  producing  a  regular  reflection 
region.  When  the  shock  reaches  a  ground  range  approximately  equal  to  the 
HOB  an  abrupt  transition  to  Mach  reflection  occurs.  This  transition  is 
responsible  for  an  airblast  environment  more  severe  than  the  surface  burst 
nuclear  case.  Qualitatively,  it  can  be  thought  of  as  a  partial  flow  stagnation 
in  the  Mach  region  that  leads  to  the  production  of  two  static  pressure  peaks. 

A  1  Kiloton  (1  K.T)  atmospheric  nuclear  explosion  at  a  HOB  of  104  feet  has 
been  simulated  using  the  two  dimensional  FAST2D  code  (Ref  .1).  Figure  1 
illustrates  the  shock  structure.  The  calculation  predicts  the  transition  of  the 
shock  from  regular  reflection  to  double  Mach  reflection.  Because  the 
spherical  waves  are  expanding  and  thus  decreasing  in  Mach  number  as  well  as  angle  of 
incidence  with  the  ground,  they  create  a  dynamic  Mach  stem  formation.  In  comparision 
to  planar  shocks  on  wedges  one  finds  them  to  be  qualitatively  alike.  The  appearance  of 
double  peaks  in  the  pressure  and  density  profiles  (versus  time  and  distance)  is 
interpreted  as  the  point  of  transition.  Other  interesting  phenomena  such  as  the 
rollup  of  the  contact  surface  generating  a  vortex  ring  and'  the  associated 
phenomenon  of  toeing  out  of  the  first  Mach  stem  can  be  observed. 

The  ability  of  the  calculation  to  accurately  predict  the  gasdynamic  effects 
both  temporally  and  spatially  is  due  in  part  to  the  shock  capturing  and  adaptive 
rezone  features  of  the  FAST2D  code.  A  minimal  number  of  very  fine  zones  was  placed 
around  the  shock  front  and  these  zones  then  moved  with  the  first  Mach  stem  to 
prevent  shock  smearing  and  distortion.  This  calculation  is  the  first  attempt  to 
model  the  nuclear  HOB  case  through  the  use  of  a  Flux-Corrected  Transport  (FCT) 

algorithm  (Ref.  2). 

Manuscript  submitted  August  24,  1981. 
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Tlie  results  of  this  calculation  agree  well  with  the  pressure-distance 
curve  generated  by  the  high  explosive  (HE)  data  of  Carpenter  (Ref.  3)  and  the 
analysis  of  Kuhl  (Kef.  4).  The  peak  overpressure  of  the  first  shock  at  .he  time 
of  transition  is  about  4300  psi.  Our  simulation  was  run  to  11.6  ms  (total 
time  with  t  =  3.76  ms),  which  corresponds  to  pressure  peaks  of  about  2000  psi. 
In  the  regular  reflection  region  the  peak  values  tend  to  be  about  20%  low  due  to 
the  clipping  of  the  FCT  algorithm  and  inaccuracies  in  the  initialization  of  the 
flow.  Reducing  the  minimum  zone  size  from  5  cm  to  1  cm  in  a  one-dimensional 
test  calculation  eliminated  this  discrepancy,  however.  In  the  two-peak  regions 
the  agreement  between  the  experimental  data  and  the  values  presented  here  is 
very  good.  The  resolution  of  the  calculation  is  adequate  for  studying  quali¬ 
tatively  the  characteristics  of  the  flow  field.  For  future  work  we  recommend 
that  the  transition  region  be  explored  with  improved  resolution. 
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II.  DESIGN  OF  PROBLEM 


The  problem  of  a  1-K.T  nuclear  detonation  at  104  ft  (31.7m)  HOB  was  chosen 

since  it  can  be  scaled  conveniently  to  various  HE  tests.  The  use  of  the  1KT 

standard  is  also  expendient;  one  could,  however,  have  used  realistic  initial 

conditions,  such  as  the  Los  Alamos  Scientific  Laboratory  RADFLO  or  Air  Force 

Weapons  Laboratory  (AFWL)  SPUTTER  calculations.  A  simple  constant  ambient 

-3  3 

atmosphere  was  used  with  a  density  of  1.22  x  10  g/cm  and  pressure  1.01  x 
6  2 

10  dyne/cm  .  To  relate  the  energy  and  mass  densities  to  the  pressure, 
a  real-air  equation  of  state  (EOS)  was  used.  This  "table-lookup"  EOS  is  derived 
from  Gilmore's  data  (Ref.  5.)  and  has  been  vectorized  for  the  TI  Advanced 
Scientific  Computer  at  NRL  (Ref.  6).  Figure  2  illustrates  the  effective  gamma 
versus  specific  energy  per  unit  mass  for  different  values  of  the  density.  The 
internal  energy  density  used  in  the  call  to  the  EOS  is  found  by  subtracting 
kinetic  energy  from  the  total  energy;  this  can  be  negative  due  to  phase  errors 
in  the  fluid  variables.  When  this  occurs,  the  value  of  the  pressure  is  reset  to 
zero. 

The  transition  from  regular  reflection  to  double  Mach  reflection  is  known  to 
occur  at  a  ground  range  approximately  equal  to  the  HOB.  Therefore,  the  size  of  the 
mesh  should  be  roughly  twice  the  HOB  in  both  directions.  The  upper  boundary 

should  be  far  enough  away  from  the  blast  front  to  be  noninterfering.  We  set  the 

3  4 

boundaries  at  5.5  x  10  cm  for  the  radial  direction  and  1.035  x  10  cm  for  the 

axial  direction.  The  fine  grid  in  the  radial  direction  contained  140  out  of 

200  total  zones,  each  5  cm  in  length.  The  largest  zones  initially  filled  the 

right  section  of  the  grid  and  were  80  cm  in  length.  A  smoothing  involving  40 

zones  was  performed  between  the  region  to  guarantee  that  the  zone  sizes  varied 

slowly.  In  the  vertical  direction  the  fine  grid  contained  75  out  of  150  total 

zones,  each  5cm  in  length.  Beyond  that  region  the  zones  increased  geometrically 

by  a  factor  of  1.112. 
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Placement  of  the  fine  grid  at  the  origin  (ground  zero  -  the  point  at  which 
first  reflection  occurs)  was  determined  to  be  optimum  for  capturing  peak  pressures 
in  the  airblast  wave  front.  Thus,  as  the  expanding  wave  moved  along  the  ground 
surface,  the  fine  grid  was  always  locked  to  it,  and  each  point  along  the  incident 
blast  front  encountered  the  same  spatial  gridding  as  it  approached  the  ground. 

By  treating  each  point  of  the  incident  front  in  the  same  manner,  we  insured  that 
the  calculation  was  internally  consistent  and  that  the  computed  transition  point 
was  accurate  to  within  the  limits  of  the  resolution.  Finally,  we  point  out  that, 
as  a  section  of  the  incident  blast  wave  propagated  within  the  fine  grid,  the 
wave  steepened.  The  size  of  the  fine  grid  was  sufficient  to  insure  that  the  incident 
wave  had  reached  the  maximum  steepeness  prior  to  intersecting  the  ground. 

The  initialization  provides  a  strong  shock  with  Mach  number  =  12.  This 
speed  and  the  need  for  restart  capability  led  to  the  choice  of  200  timesteps 
as  an  interval  for  the  spatial  display  ("snapshots") .  The  dump  interval  that 
resulted  was  ~At  =  0.3  milliseconds.  These  dumps  were  stored  on  magnetic  tape 
and  postprocessed. 

Additional  diagnostics  were  implemented  in  the  calculation.  Stations  were 
created  to  gain  information  from  fixed  spatial  positions  within  the  calculational 
grid.  These  25  physical  variable  sensors  were  placed  along  the  ground  and  stored 
values  of  the  energy  and  mass  densities  and  velocity  for  every  timestep.  From 
this  information  one  can  construct  static  and  dynamic  pressure  curves. 

III.  COMPUTATIONAL  DETAILS 

The  evolution  of  the  nuclear  HOB  flow  field  was  modeled  numerically  with  the 
FCT  code  FAST2D  (Ref.  7).  FCT  yields  accurate  and  well-resolved  descriptions 
of  shock  wave  propagation  without  the  necessity  of  a  priori  knowledge  of  the 
essential  gasdynamic  discontinuities  in  the  problem.  Additionally,  the  code  has  a 
general  adaptive  regridding  capability  which  permits  fine  zones  to  be 
concentrated  in  the  region  of  greatest  physical  interest  while  the 
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remainder  of  the  system  is  covered  with  coarse  zones.  Figure  3  depicts  the  grid 
setup  initially  and  at  transition  to  the  double  Mach  stem  structure.  The  rezone 
algorithm  is  programmed  to  track  the  Mach  stem  with  the  fine  grid. 

The  transport  algorithm  used  a  low-phase-error  phoenical  FCT  algorithm  in  a 
model  called  JPBFCT,  an  advanced  version  of  the  ETBFCT  algorithm  described  in  Ref  1. 
The  linear  part  of  this  algorithm  is  fourth-order  accurate  spatially  in  advection 
problems  with  a  given  constant  velocity  and  has  a  (nonlinear)  flux-corrected 
antidiffusion  needed  to  model  shocks  correctly.  Finally,  the  transport  subroutine 
is  written  in  sliding-rezone  form,  which  means  that  the  mesh  at  the  beginning  and  the 
end  of  the  timestep  need  not  be  the  same.  Since  the  algorithm  is  one-dimensional, 
timestep  splitting  is  employed  to  solve  the  2-D  problem. 

The  fluid  transport  routine  JPBFCT  is  fully  vectorized  and  requires  about 
2  us  per  meshpoint-cycle .  This  time  would  have  been  still  less  if  a  vectorized 
fully  two-dimensional  routine  had  been  used,  since  the  1-D  loops  are  too  short  to 
permit  full  advantage  to  be  taken  of  the  vector  capabilities  of  the  NRL  ASC.  The 
table  lookup  in  the  EOS  was  also  fully  vectorized,  so  that  pressure  calculations 
required  about  20%  of  the  time  needed  for  the  hydrodynamics.  These  two  items 

took  up  nearly  all  of  the  running  time  in  the  blast  calculation  itself.  The  cost 

of  initialization  was  negligible,  but  the  diagnositics  cost  up  to  30%  as  much  as 
the  hydrodynamics,  depending  on  how  many  of  the  various  possible  quantities  were 
actually  plotted.  This  latter  number  would  be  greatly  reduced  if  the  plot  routines 
were  fully  vectorized. 

A  version  of  the  AFWL  1  KT  standard  (Ref  8)  was  used  to  initialize  the  energy, 

density  and  velocity  (flow  field)  at  3,76  milliseconds.  The  corresponding  shock 

4 

radius  was  103.9  ft  (31.69  m)  peak  overpressure  of  1645  psi  (1.134  x  10  K  Pa). 
Because  some  areas  of  the  grid  were  very  coarse,  interpolation  onto  the  grid  was 

performed.  After  the  1  KT  flow  was  laid  down  inside  a  radius  of  104  feet  (31.7  m) 

the  fine-zoned  grid  was  activated  to  follow  the  peak  pressure  as  it  moved  along 
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Fig.  3  -  Adaptive  gridding.  The  grid  at  initialization  and  at 
transition  point  (lines  are  drawn  for  every  other  zone,  lines 
in  fine-zone  region  are  indistinguishable). 
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the  ground  surface,  modelled  as  a  perfectly  reflecting  boundary.  This  region 
comprised  140  zones,  and  a  switch  was  set  to  keep  40  of  these  zones  ahead  of 
the  reflection  point.  Permeable  boundary  conditions  were  used  on  the  top  and 
right  edges  of  the  mesh;  i.e.,  density,  pressure  and  velocity  were  set  equal 
to  ambient  preshock  conditions.  Reflecting  conditions  were  applied  to  the 
left  and  bottom. 

The  timestep  was  recalculated  at  every  cycle  according  to  the  Courant 
condition 


( Ax  t , Ay . ) 

dt  =  0.5  min  (V+fvT') 

ij  ij 


(1) 


where  c  is  the  speed  of  sound  and  jvj  is  the  modulus  of  the  flow  speed.  This 
could  have  been  relaxed  somewhat  by  allowing  violation  of  the  local  Courant 
limit  at  points  ("hot  spots")  far  from  the  region  of  chief  physical  interest. 

The  total  elapsed  time  in  the  2-D  calculation,  7.6  ms,  required  5600  cycles. 

Three  types  of  diagnostics  were  employed,  all  in  the  form  of  plots  made 
by  post  processing  a  dump  tape.  The  first  type  of  diagnostic  consisted  of  CRT 
contour  plots  of  density  and  static  pressure,  and  arrows  indicating  the  magnitude 
and  direction  of  the  velocity  field,  obtained  at  the  dump  intervals  (every  200 
cycles).  The  second  type  was  pressure-range  curves  at  z=0,  obtained  by  finding 
the  pressure  peak(s)  along  the  ground  at  each  dump  interval  and  hand-plotting 
them  on  the  same  graph.  The  third  type  consisted  of  pressure  histories  at  a  series 
of  24  stations,  obtained  by  saving  the  energy  and  mass  densities  and  the  veloci¬ 
ties  at  every  cycle. 

IV.  RESULTS  AND  PHENOMENOLOGY 


This  calculation  has  been  done  to  understand  the  violent  effects  of  1  KT 
of  energy  being  released  in  the  atmosphere  at  a  HOB  of  104  ft  (31.7m).  A  strong 
spherical  shock  is  created  in  the  surrounding  air,  and  reflects  from  the  ground. 
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The  outward-traveling  airblast  is  then  conposed  of  two  parts:  one  reflected 
upward  approximately  normal  to  the  ground,  and  the  original  spherical  blast. 

The  peak  pressure  is  coincident  with  the  intersection  of  the  two  waves.  This 
intersection  continues  to  move  outward  until  the  angle  of  the  spherical  shock 
with  respect  to  the  ground  reaches  a  critical  value  and  the  transition  to  a 
double  Mach  stem  occurs.  As  shown  by  Ben-Dor  and  Glass  (Ref.  9),  this  angle 
depends  upon  the  incident  strength  of  the  shock  (Fig.  4).  Shocks  with  Mach 
numbers  greater  than  10  are  not  shown.  Initially,  the  Mach  number  for  the  HOB 
simulation  is  well  above  10.  The  Mach  number  at  transition  is  approximately 
11  and  the  angle  is  less  than  50°.  From  Fig.  4  the  corresponding  region  is 
double  Mach  reflection. 

Figure  1  has  been  labeled  with  the  notation  of  Ben-Dor  and  Glass  (Ref.  9). 

It  should  be  noted  that  what  is  generally  regarded  as  the  second  Mach  stem 
is  in  fact  the  second  reflected  wave,  which  is  part  of  the  second  Mach  struc¬ 
ture.  To  be  consistent,  one  must  label  the  second  Mach  stem  as  at  the 

indicated  location.  The  definition  used  is  the  state  of  the  fluid  one 
obtains  by  passing  through  one  shock  wave  (M^)  or  two  shock  waves  (R  and  R^). 

The  first  reflected  wave  R  becomes  the  incident  wave  for  the  second  Mach 
structure.  Density  contours  are  shown  in  Fig.  5  for  an  planar  shock  on 
wedge  with  a  Mach  number  of  7  and  an  angle  of  50°.  The  complimentary 
figure  illustrates  the  proper  labeling  of  the  multiple  waves.  Comparison 
of  Fig.  5  and  the  HOB  simulation  (Fig.  6)  shows  that  corresponding  waves 
can  be  identified.  Differences  between  the  planar  shock  on  wedge  and  the 
HOB  can  be  explained  in  terms  of  the  unsteady  nature  of  the  HOB  case  (a 
spherically  expanding  wave  that  continuously  decreases  in  Mach  number  and 
angle.)  Although  the  term  irregular  Mach  reflection  has  been  used  to 
describe  the  complex  shock  structure  that  evolves  from  HOB  events,  we  believe  it 
to  be  very  regular  and  explainable  as  a  double  Mach  reflection  that  evolves  as  a 


function  of  time. 


WEDGE  ANGLE  (DEGREES) 


90 


MACH  NO.  - - 

Fig.  4  -  Types  of  shock  reflection:  RR,  CMR,  SMR,  and  DMR  denote 
regular  reflection  and  single,  complex  and  double  mach  reflection, 
respectively.  The  D  and  A  refer  to  attached  and  detached  shocks. 
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with  wave  identification 


The  HOB  numerical  simulation  begins  just  before  the  shock  first  reflects 
from  the  ground.  As  a  summary  of  how  the  flow  field  then  develops,  we  pre¬ 


sent  snap-shots  at  the  important  stages .  (A  more  complete  display  is  presented 
in  Appendix  A) .  Figure  6a  indicates  the  pressure  and  density  contours 
and  velocity  vectors  at  t  =  3.18  ms.  In  Fig.  6b  the  reflected  shock  is  shown 
moving  upward,  the  outward  flow  begins  to  stagnate  at  the  ground  (transition). 

Fig.  6c,  at  t  =  5.99  ms,  shows  an  enlargement  of  the  shock  front,  and  the  development 
of  the  Mach  stem,  slip  surface  and  second  Mach  stem.  The  angle  of  the  shock 
with  respect  to  the  ground  is  increasing  with  time,  so  that  the  effective  wedge 
angle  is  decreasing.  From  the  work  of  Ben-Dor  and  Glass  one  expects  a  transition 
to  double  Mach  stem  to  occur  at  approximately  45°.  The  angle  in  Fig.  6b  is  about 
45°  and  the  shock  front  has  entered  the  transition  phase.  Figure  6d  shows  the 
fully  developed  shock  structure  at  7.79  ms.  Toeing  out  of  the  first  Mach  stem 
can  be  also  seen  in  the  contours  of  Fig.  6d  and  occurs  as  the  fluid  rolls 
forward  where  the  slip  line  would  otherwise  intersect  the  ground.  The  velocity 
field  in  Fig.  6d  also  shows  this  detail. 

Note  the  reflected  shock  properties  (that  part  of  the  structure  that  contains 

the  second  Mach  stem  M^).  The  reflected  shock  propagates  rapidly  through  the  high 

temperature  fireball,  due  to  the  high  local  sound  speed.  The  shape  of  this  reflected 

wave  is  a  primary  difference  between  the  HOB  case  and  the  planar  wave  on  wedge  case. 

The  other  major  difference,  of  course,  is  the  spherically  expanding  blast  wave  which 

-2 

decreases  in  strength  approximately  as  r 
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An  effective  way  to  quantitatively  evaluate  the  calculation  and  observe 
in  detail  the  transition  to  a  Mach  stem  regime  can  be  seen  by  examining  the 
station  data.  The  station  sensors  were  placed  in  the  bottom  row  of  the  cal¬ 
culation  100  to  200  cm  apart.  In  Table  I  the  maximum  pressure  recorded  for 
each  station  along  with  the  location  can  be  found. 

Besides  giving  a  reliable  value  for  the  peak  pressure  to  be  used  for 
constructing  the  pressure-distance  curve,  these  data  allow  one  to  see  effects 
fixed  in  space  but  varying  in  time.  Figure  7  is  a  superposition  of  the  pressure 
profiles  from  stations  15  to  24  (the  transition  region).  Noteworthy  is  the 
profile  from  station  17,  which  is  the  first  station  to  record  a  second  peak  on 
the  back  side.  At  station  19,  200  cm  further  away  from  ground  zero,  the  second 
peak  is  almost  equal  to  the  first.  The  visible  transition  (as  seen 
in  Fig.  6b,  occurs  at  a  ground  range  between  3100  cm  and  3300  cm  (stations 
19-21)  revealing  a  dominant  second  peak  and  a  "first  peak"  (i.e.,  first 
seen  by  the  sensor)  that  is  about  half  the  magnitude  of  the  second.  The 
second  peak  does  not  exhibit  a  sharp  almost  discontinuous  rise  and  then  a 
rapid  but  slower  decrease  along  the  back  side.  Instead,  it  has  the  appearance 
of  a  density  compression.  This  behavior  has  dramatic  consequences  for  mi¬ 
litary  planners  because  the  pressure-distance  curve  is  modified  and  the 
dynamic  pressure  is  enhanced. 

The  analogous  profiles  for  dynamic  pressure  are  presented  in  Fig.  8. 

Again  data  from  stations  15  to  24  is  utilized.  The  development  of  the 
second  peak  and  its  correlation  with  the  Mach  stem  formation  can  be  observed. 
There  is,  in  addition,  a  noticeable  increase  in  the  first  peak  values  (station 
15  to  the  maximum  at  station  18).  After  the  structure  becomes  visibly 
resolved  (station  20  and  beyond)  the  second  peak  resembles  a  rounded  profile 
suggesting  the  formation  of  a  stagnation  region  behind  the  first  peak  (Mach 
stem) . 
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Table  1 


Station  No. 

Location  (cm) 

Time  (sec) 

Pres  (dynes/ 

1 

2.0000E  02 

2.25E-0A 

8. HE  08 

2 

A.0000E  02 

2.80E-0A 

7.92E  08 

3 

8.0000E  02 

5.28E-0A 

7.17E  08 

A 

1.0000E  03 

7.23E-0A 

6.73E  08 

5 

1.2000E  03 

9.5AE-0A 

6.2AE  08 

6 

1.A000E  03 

1.23E-03 

5.65E  08 

7 

1.6000E  03 

1.5AE-03 

5.21E  08 

8 

1.8000E  03 

1.91E-03 

A.70E  08 

9 

2.0000E  03 

2.3AE-03 

A.5AE  08 

10 

2.2000E  03 

2.81E-03 

A.1AE  08 

11 

2.3000E  03 

3.07E-03 

A.03E  08 

12 

2.A000E  03 

3.35E-03 

3.92E  08 

13 

2.5000E  03 

3.62E-03 

3.82E  08 

1A 

2.6000E  03 

3.88E-03 

3.73E  08 

15 

2.7000E  03 

A.19E-03 

3.37E  08 

16 

2.8000E  03 

A.A8E-03 

3.33E  08 

17 

2.9000E  03 

A.79E-03 

3.05E  08 

18 

3.0000E  03 

5.11E-03 

3.01E  08 

19 

3.1000E  03 

5.A1E-03 

2.33E  08 

20 

3.2000E  03 

6.06E-03 

1.96E  08 

21 

3.3000E  03 

6.A9E-03 

1.79E  08 

22 

3.A000E  03 

6.82E-03 

1.87E  08 

23 

3.5000E  03 

7.28E-03 

1.85E  08 

2A 

3.6000E  03 

7  .73E-03 

1.69E  08 

Finally  we  consider  the  pressure-distance  relation  for  the  HOB  case. 

In  Fig.  9  we  compare  the  results  of  the  numerical  simulation  with  the  data 
of  Carpenter  and  with  empirical  analysis.  Carpenter's  data  are  based  upon 
careful  HOB  experiments  with  8  lb  PBX9404  spheres.  The  empirical  analysis 
was  based  on  a  1  KT  nuclear  free  air  curve  and  HOB  construction  factors. 

The  calculated  values  in  the  regular  reflection  regime  are  20%  low,  which 
may  be  attributed  to  a  combination  of  FCT  clipping,  the  resolution  of  the  grid, 
and  inaccuracies  in  the  initialization  of  the  flow  field.  During  and  after  Mach 
reflection,  the  peaks  remain  low  until  the  Mach  stem  structure  has  grown  large 
enough  to  be  resolved  on  the  mesh.  By  the  time  it  occupies  a  region  of  15  cells 
high  and  35  cells  wide,  the  peak  pressures  are  in  good  agreement  with  the  HE 
data  and  the  empirical  analysis. 

Other  attempts  to  model  the  transition  region  have  been  made.  Needham  and 
Booen  (Ref.  10)  present  results  of  a  1100  lb  pentolite  sphere  at  15  feet  HOB.  The 
general  phenomena  of  the  flow  field  can  be  seen  from  their  simulation.  When  a 
pressure  distance  curve  is  constructed  from  this  calculation,  one  finds  that  in  the 
regular  reflection  region  their  results  are  15%  to  30%  high  relative  to  theory. 
After  transition  to  double  Mach  reflection  the  first  peaks  are  20%  low  while 
the  second  peaks  are  40%  low  (Ref.  4). 

V.  SUMMARY 

The  airblast  from  a  1KT  nuclear  event  at  104  ft  HOB  has  been  numerically 
simulated  with  the  FAST2D  computer  code.  The  results  give  insight  to  the 
formation  and  subsequent  evolution  of  the  Mach  stem,  the  triple  point,  and 
the  contact  discontinuity.  The  transition  from  regular  reflection  to  double 
Mach  reflection  is  predicted.  We  suggest  that  the  first  signal  for 
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transition  is  the  appearance  of  a  second  peak  behind  the  shock  front  due  to 
stagnation  in  the  flow.  Comparison  with  the  pressure-distance  curves  of 
Carpenter  and  Kuhl  indicates  agreement  within  20%.  Both  first  and  second 
peaks  are  predicted  with  similar  accuracy. 
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APPENDIX  A.  DETAILED  TIME  HISTORY  OF  CALCULATION 


The  following  figures  comprise  a  temporal  history  of  the  numerical 
simulation.  Each  page  contains  pressure  contours,  velocity  vectors, 
density  contours,  and  the  corresponding  grid  for  a  particular  time. 

The  series  begins  at  t  =0  (tT  =  3.76  ms)  and  continues  to  t  =  8.28  ms. 
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OICY  ATTN  AtJl-LSE 


rp CAR TMFNT  PF  THE  M  p  FCPCF 


ASSISTANT  CHIEF  rr  STAFF 

intelligence 

DEPARTMENT  OF  The  AIc  cC*Cr 
WASHINGTON,  DC 

0 1C  Y  ATTN  IN  cm  4  AS?? 


ASSISTANT  CHIEF  CF  STAFF 
STUDIES  C  ANAL YS  FS 
PEoA0TMFNT  of  THF  A IO  FPOCE 
WASHINGTON,  DC  ?0??0 

0  ICY  ATTN  AF/SA^I  <  T f  CH  Hr}} 


ASSISTANT  SECFFTApy  OF  THE  AF 
PF SEARCH,  DEVELOPMENT  p  LOGISTICS 
DEPARTMENT  OF  THF  A  I  0  PC  PC  F 
WASHINGTON,  DC  2???'' 

DICY  ATTN  SAP Alp /CE°  FOR  $TRAT  p  SPACE  SYS 


BALLISTIC  MISSILF  OFFICF/on 
air  FOPCF  SYSTEMS  CC"*ANC 
NORTON  AFQ,  CA  9?4CS 
( M  INiJTE  MAN  ) 

OICY  ATTN  MNNXH  G  KALANSKY 

01CY  ATTN  MNVVH  V  CFLVFCC.HIO 

OICY  ATTN  WNN  W  CFADTREE 

OICY  ATTN  MNNXH  D  GAGF 

OICY  ATTN  M.NNY 


DEPUTY  CHIEF  hp  STACF 
oESFARCH,  CEVFICPVEhT,  p 
DFpARTMcmt  or  i>r  4  Tp  crpcr 

WASHINGTON,  DC  ?0??G 

OICY  ATTN  AFcnui  \  ALFXANDROW 
OICY  ATTN  A*"p  pp'l 
OICY  ATTN  AFPOQI 
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DFPART'-'FNT  CF  THE  AlP  FCPCE 


DEPUTY  CHIEF  OF  STAFF 
LOGISTICS  C  ENGI  NFr;R  ING 
DEPAP  TMENT  OF  THE  AI°  TCRCr 
WASHINGTON,  DC  20  3  3C 
OICY  ATTN  LEFF 


COMMANDER 

FOREIGN  TECHNOLOGY  CIVISIQN,  AFSC 
WR I GH T— P ATT EP S ON  AFP,  PP  45 ATT 
OICY  ATTN  Ml  I S  L  I  PP ARY 


COMMANDER 

P?'>F  AIR  DEVELOPMENT  rc^T^^t  AFSC 
GRTFFISS  AFB ,  Nv  13441 

(CESIPFS  NC  CM’-'OT) 

OICY  ATTN  TSLO 


STRATEGIC  AIR  CQMVANC 
DFPAPTMENT  OF  THF  A  IP  PPPCE 
C.PrUTT  AFB,  NB  6R  1  1  * 

OICY  ATTN  NPI-STINFO  LI«°ARv 

OICY  ATTN  XPFS 

OICY  ATTN  TNT  J  NO'IM'.EY 


VELA  SE I S^OLDGI  0.  Al  CFf'TrP 
312  MONtg^vcry  S TRF  F  t 
ALEXAMDIRA,  V A  2? 314 
OICY  ATTN  G  "i (  o  ICF 


HEPAOTMfMT  or  ENERGY/DOF  CONTRACTORS 


DEPARTMENT  or  ENERGY 

AL  3UCUFFF'  if  pr>cnATICAS  C^Ifr 
P  0  BOY  5A00 
AL SUOl'C^O1  |C  *  VM  F7115 
OICY  ATTN  C T T n 


PEP  AR  TmFNT  of  r\>cpr,Y 
WASHINGTON,  HC  70545 
OICY  ATTN  C’-1 A  /PD£  7 


DEPARTMENT  OF  EN FP GY 
NEVADA  nP^PATIrNS  OFFICE 
P  0  3°X  l A 100 
LAS  VEGAS,  NV  39114 

OICY  ATTN  f,A  I  L  £  RECORDS  FOP  TECHNICAL  LIPPARY 


LAWPEMCF  l  IVE°nppf  NAT  I C>  Al  t AB 
P  0  BOX  POR 
LIVERMORE,  CA 

OICY  ATTN  !-nC  P  DC N G 

OICY  ATT N  L -205  J  FEA°ST  (CLASS  L~?0?) 

OICY  ATTN  L  PD  ^  rrPRJS  (CLASS  L-594) 

OICY  ATTN  L-7  J  KAHN 

OICY  ATTN  D  01  CNN 

OICY  ATTN  L  A  3  7  R  SCHCCK 

OICY  ATTN  TFCFNTCAI  IN^O  ncn.  Ll^ARY 

OICY  ATTN  l -7  CD  T  Bl TKDV I CH 


LOS  ALAMOS  NATIONAL  SriFNTJPTC  LAB 
MA  I  L  ST  AT  I  ON  r-000 
P  O  BOX  1463 
LOS  ALAMOS,  NM  P  7  5  4  5 

(CLASS! FIFO  rNl Y  Tf  MAIL  STATION  5000) 
OICY  ATTN  P  WMTTA*FR 
-Q-K Y  -ATTN-P-  C--FOV  Cl 4- 
0  ICY  ATTN  C  KEIL<=R 
OICY  ATTN  PTANf?n  d  ■  ^  T 
OICY  ATTN  vs  it'-  (CLASS  pf°°RTS'LIB) 

0  i  Lrl  7TTO  P=  J,_  *j  .3.  S' 
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DEPARTMENT  OF  ENERGY/DOF  CONTRACTORS 


LOVELACE  PIOMFOICAL  F 

FNVIPONMFNTOL  RSCH  P'SHU1^  INC. 
P  0  Rpy  3890 
ALBUQUFPQUF,  N*  FT)  15 

OICY  ATTN  »  JCNCS  (UNCI.  CMY) 


OAK  PTOCF  NAT  10*’ Al  L  A  00° A  TOP  V 

NICLFAP  DIVISION 

X~in  LAD  PE CNR OS  DIVISION 

P  0  ROX  X 

CAY  Pinr.F,  TM 

OICY  ATTN  CIVIL  1 EF  PFS  PPOJ 
OICY  ATTN  CErT°A(  "SCH  LIBRARY 


SANDIA  LABORATORIES 
LIVER  YORE  L  AP,  QR  A  tdo  y 
p  0  RnX  Dfcq 

LI  vcp.Mnpp  ,  CA  °AC5D 


SANDIA  NATIONAL  I  AR 
P  0  COX  5P0Q 
ALBUQUERQUE ,  NM  87139 


OICY 

(ALL 

ATTN 

CLASS  ATTN  s^c 

A  CHAR  AN 

OICY 

ATTN 

t  HILL 

OICY 

ATTN 

OPG  I? 50  Vj  P^OwL 

OICY 

ATTN 

A  CI-A3  I A 

OICY 

ATTN 

W  POHEPTY 

OICY 

ATTN 

314! 

OICY 

ATTN 

L  V OPT MAN 

Ol  v ^ 

/yr*J 

■j?  st’i^r  / 
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nFPAPTVFNr  CF  DEFENSE  CCNTPACTQBS 


AC'IPEX  rr^P. 

435  CLYD r  AVFvijf 
MOUNT  A IV  VIEW,  CA  0  404? 
OICY  ATTN  C  vrLF 


AEROSPACE  CO°°. 

P  C  BOX  9?P57 

LOS  ANGFL  ec  ,  CA  9000^ 

OICY  ATTN  H  ”IPG19 

OICY  ATTN  TECHNICAL  I NFCP NATION  SERVICES 


AGBABIAN  ASSOCIATE 
250  N  NASH  ST ° EE T 
EL  SFGMNOn  ,  CA  90?45 
OICY  ATTN  m  AGBABIAN. 


ANALYTIC  SF°VICES,  INC. 

400  APMY-NAVY  ro  j\/c 

Arlington,  va  22202 

OICY  ATTN  o  MESSFLBACHEP 


APPLIED  0 E S EApC H  ASSCCIATES,  INC 
2601  WYOMING  BLvn  NF  S"ITE  H-l 
ALBUOUEPOME,  N"'  P71I? 

OICY  ATTN  J  p  r  A  T  T  CN 
OICY  ATTN  N  HI  GO  I  VS 


A°2  L I EO  THEORY,  INC. 

1010  WFSTwpnn  plvt 
LOS  ANGELES,  CA  90.  n?  4 

(2  CYS  IE  l NC L ASS  op  \  CY  IF  CLASS) 
OICY  ATTN  J  Tout  rr 


pcoAPT^cvr  nr.  OFFENSE  CONTRACTORS 


A»Tcp  ASSOCIATES,  IAC. 
>6046  F9FN  LADING  P^AO 
HAYWAPn,  CA  94545 
OICY  ATTN  9  3  II  L 


ASTRON  PESAPCH  £  ENG  INF EF  ING 
1091  010  MIOOLFF  ITI.C  WAY  *15 
VOJNTAIN  VIEW,  CA  9  494  3 
OICY  ATTN  J  FONT T NO  TON 


AVCO  RESEARCH  £  SYSTEMS  CRptip 
?0l  LOWFLl  STREE T 
WILMINGTON,  MA  91B97 

OICY  ATTN  t.IHRAPY  AF30 


BOM  COPP. 

7915  JONES  BRANCH  OPIVE 
MCL  FAN ,  VA  ??I9? 

OICY  ATTN  4  IAVAGMMC 
OICY  ATTN  T  NEIGHPEP S 
OICY  ATTN  COR Rf? RATE  I  !  A  P v 


BOM  rnpo, 

P  0  BOX  9?7 4 
ALBUOUFPO"E,  NM  87119 
OICY  ATTN  R  HENSLFV 


BOEING  CO. 

P  0  BOX  3707 
SEATTIF,  W A  0*174 
OICY  ATTN  S  STRACK 
OICY  ATTN  AFRCSRACF  LIBRARY 
OICY  ATTN  M/S  4?/>7  9  C A° l  SON 
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DEPARTMENT  OP  DEFENSE  CCNTR ACTORS 


CALIFORNIA  RESEARCH  5  TECHNOLOGY,  INC. 
6269  VA R I c L  AVENL'E 
WOODLAND  HILLS,  CA  01367 
01CY  ATTN  L I RPA°  v 
01CY  ATTN  K  KREYENHAGEN 
0  ICY  ATTN  M  R0SFN61 ATT 


CALIFORNIA  RESEARCH  C.  TECHNOLOGY,  INC. 
4049  FIRST  STPFFT 
L  IVFP  MORF,  CA  C4ESC 
OICY  ATTN  D  PP°HAL 


CAL  SDAN  COR  P. 

P  0  ROX  4°0 
9’IFFALC,  NY  1422  0 
OICY  ATTN  LIPP^Y 


OEMVFP,  UNIVERSITY  OF 
COLORADO  SEMINARY 
OFNVFR  RFSFAPCH  INSTITUTE 
P  0  BOX  10127 
DENVFR,  CO  80210 

I  ONLY  1  COPY  OF  CLASS  RPTS) 
OICY  ATTN  SEC  PFFICFR  FOR  J  WISOTSKI 


EGCG  WASH.  ANAtYTICAl  SVCS  CTR,  INC. 
P  0  PPX  1021 p 
ALR'JOUEROUE,  Nv  87  114 
OICY  ATTN  LIBRARY 
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DEPARTMENT  OP  OEFENST  CCNTP  ACTORS 


ERIC  H.  WANG 

CIVIL  ENGINEER  INC  -SCF  FAC 
UNIVERSITY  OF  '-Fyicr 

UNIVERSITY  ST AT  T ON 
P  n  BOX  25 

ALRUCUEPQUE,  NM  R 7  1  “*  l 
OICY  ATTN  J  l  AY 0 

o ic y  attn  o  trmc 

OICY  ATTN  M  3 A' '*•' 

OICY  ATTN  J  KCVARNA 


GARD,  INC. 

7449  N  NATCHEZ  AVEVt  E 
NILES,  IL  60640 

OICY  ATTN  G  YC!nFAP6T  I'lNCl  ONLY) 


GENERAL  ElFCTGIC  CO. 

SPACE  OT  V I S InN 
VALLEY  FORGE  SPACE  CEf'TEP 
P  0  enx  R555 
PHILADELPHIA,  94  19101 

OICY  ATTN  M  30PTNFP 


GENFRAL  RESEARCH  CHRP. 
SANTA  PARR A PA  OTVISTCN 
P  0  30X  6770 
SANTA  PAP.BARA,  CA  92111 
OICY  ATTN  Tjr 


M-TFCH  LABS,  INC. 

P  n  BOX  1686 
SANTA  MONICA,  CA  90406 
OICY  ATTN  B  HAPTENPALV 


DERARTMFNT  OF  OFFFN'SF  CENT  RAC  TOPS 


HORIZONS  TECHN^ICGY,  TNC. 
7830  CL  A  IP  E  MONT  N'ESA  BLVD 
SAN'  01  EGO,  CA  "Pill 
OICY  ATTN  R  KPMGFP 


1 1  T  RESEARCH  TNSTTTLTC 
10  W  ?5TH  STPRCT 
CHICAGO,  TL  60616 
0  IT Y  ATTN  P  1  FLCF 
OICY  ATTN  ,J  JOHNSON 
OICY  ATTN  OOCI'YONTS  IIP&APY 


INFORMATION  SCIFNCF,  INC. 
123  W  OAORF  ST°EFT 
SANTA  BARBARA,  CA  93105 
OICY  ATTN  W  Olio?  i/jk 


INSTITUTE  EOF  DEFENSE  ANALYSTS 
400  ARMY-NAVY  OP  IVF 
ARLINGTON,  V A  2220? 

OICY  ATTN  CLASSIFIED  LIBRARY 


J  D  HALTIWANGER  CCNSLLT  ENG  SVCS 
RM  106A  CIVIL  ENGINEERING  Bl  CG 
208  N  ROM  INF  STREET 
"RRANA,  XL  A 1 RQ 1 
OICY  ATTN  W  HAH 


J.  H.  WIGGINS  CC.,  INC. 

1650  S  PACIFIC  TCAST  E IGF WAY 
REDONDO  BEACH,  CA  SC?77 
OICY  ATTN  J  COIL  INS 


59 


nFPA9TVF\'T  nr  DEFENSE  CONTRACTORS 


KAMAN  AVI  DYNE 
83  SECOND  AVFNHF 
NORTHWEST  INDUSTRIAL  PARK 
BURLINGTON,  M  A  01 -*0? 


OICY 

ATTN 

P 

R  l.  E  T  E  M  K 

OICY 

ATTN 

LIBRARY 

OICY 

ATTN 

M 

HPPB  E 

OICY 

ATTN 

C 

C  P  I  S  C I  f  N  E 

KAMAN  SCI CNCFS  C.rpp. 

P  0  POX  7 A 63 

COLOR  AOO  SPPINOS,  CC  8C933 
OICV  ATTN  0  SACHS 
01CY  ATTN  F  SHE L TCN 
OICY  ATTN  LIBRARY 


KAMAN  TEMPO 

816  STATE  STREET  (R  0  CPAUEP  00) 
SANTA  BAPPAPA,  CA  93JO? 

OICY  ATTN  DAS  I  AC 


LOCKHEED  missiles  S  SPACE  CC.,  INC. 
P  0  nPX  504 
SUNNYVALE,  CA  94086 
OICY  ATTN  J  W  F I  S  N  F  P 
OICY  ATTN  TIC-LiePARY 


MARTIN  M4DIFTTA  Cr)RR. 
P  0  BOX  5937 
OP.LANOO,  FL  33*55 
OICY  ATTN  G  FCTIFC 


MARTIN  MAPiFTTfi  fopp. 
o  0  BOX  17R 
OENVEP,  CO  °0?o ] 

OICY  ATTN  G  F°Fvr(i 
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CCpAMTMPnT  OF  DEFENSE  f.CNTRACTORS 


VCOCVNCU  DCUGLAS  CCPp. 
5101  DniSA  AVFNt :  E 
HUNT  I  NGTrN  BEACH,  r  A  c?647 
01CY  ATTM  H  HEROMA.v 
0 1C Y  ATTN  P  H  AL  °P  IK 
0  1C  Y  ATTM  0  OF  AN 


MCOOVNcLL  nn'IGLAS  crco. 
1955  1  AKFW'nr)  nn;!(_Pv*P" 
f.ONO  8CACH,  CA  <70  *47 
01CY  ATTN  m  PCI T  pp 


MFRP1TT  CASES,  INC. 

P  0  BOX  1 ? 0 6 
REDLANDS,  CA  0?371 
01CY  ATTN  J  MPT  ITT 
0 If  Y  ATTN  L  I  TP  AR  Y 


ME  T  FPROl.  CGY  RESEARCH,  J\C. 
464  w  wnnnRu°Y  pr ao 
ALT  A  OF*1  A,  CA  9100J 
01CY  A T J N  W  Gnrr\ 


MISSION  RESEARCH  CORO. 

P  0  OP  AWFC  7] 9 

SANTA  BARBARA,  CA  93  10? 

(All,  CLASS:  ajtn:  s  rc  nFC  FOP) 
OICY  ATTN  C  LCNG'MP  r 
OICY  ATTN  G  MfC  A  c  T  C  c- 


PACIFIC-SIERRA  RESEAOCH  Copn. 
1456  CL0VERFIF1.  0  Bt.VC 
SANTA  MCNICA,  CA  904C4 
OICY  ATTN  H  3D0DE 


-*  ^  i  rSl^ri  ,  ii  1 1 
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ppPAojMFMT  OF  OFFENSE  CCNTP ACTORS 


PACIF TC-S1ERRA  PfSEARCH  COpP. 
WASHINGTON  OPERATIONS 
140  1  WILSON  3L''0 

suite  uoo 

ARLINGTON,  VA  ??2  00 
Olfv  ATTN  0  r.ORMLEY 


PACIFICA  TECHNOLOGY 
P  n  POX  143 
OEL  OAR,  CA  Q?^14 
0 1C  Y  ATTN  R  B JORX 
0  ICY  ATTN  G  KENT 
OICY  ATTN  TEOH  L.  IPRARY 


PATFL  FNTrRPR I S  rS ,  INC. 
P  0  BOX  0531 
HUNTSVILLE,  AL  25R1C 
OICY  ATTN  M  °ATEL 


PHYSICS  IMTEBMAT  IONAl  CO. 

2700  v-ppcFO  STc'rrT 
SAN  I.EANOPn,  CA  c4577 
OICY  ATTN  L  p  EHR  M  ANN 
OICY  AT  TO  TECHNICAL  LIBRARY 
OICY  ATTN  «=  ’OORE 
OICY  ATTN  .1  T  h‘r  S  E*- 
OICY  ATTN  F  CAilFR 


p  £  r,  ASSOCIATE 
P  0  POX  *>605 

MARINA  CEL  PTv,  CA  OC201 


OICY 

ATTN 

p  r>np  t 

OICY 

ATTN 

A  Kl'ut 

OICY 

ATTN 

j  LF'  TS 

OICY 

ATTN 

W  0  T  “>  H  T 

OICY 

ATTN 

J  tAc>r,FVTrP 

OICY 

ATTN 

Trrp  T  C  AL  INFORMATICS  CENTER 
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CPo^pTwcnjT  CP  OFFENSE  CONTRACTORS 


PANG  COPP. 

1700  mat 'I  STRUTT 
SANTA  mpmICA,  r a  c  0  4  C  6 
OICY  ATTN  C  '*rv 


SCIENCE  APt>LIC'T  T?''S,  INC 
RAOIATION  IMST^'JVFMAT  r c ^  0 1  v 
4615  HAWKINS,  ME 
ALB'JQ'JFPQUE ,  NM  ci?irc, 

OICY  ATTN  J  niSWCN 


SCIENCE  APPLICATIONS,  INC. 

P  0  BOX  2351 
LA  JOLLA,  CA  520^o 
OICY  ATTN  H  -III  SC  N 
OICY  ATTN  TECHNICAL  L  TB0  A»Y 
OICY  ATTN  <?  SCMLALT 


SC  I FNCE  AP°L  IC  AT  ION  S  ,  INC. 
101  CONTINENTAL  Pl'/C 
FL  SFGHNOO,  CA  OC74  5 
OICY  ATTN  0  HOVE 


SCIENCE  A00LIC4TnNS,  INC. 
2450  WASHINGTON  AVENUE 
SAN  LEANDRO,  CA  04577 
OICY  ATTN  D  BERNSTEIN 
OICY  ATTN  n  MAXWELL 


SCIENCE 

APPLICATIONS,  INC. 

P  0  BOX 

130? 

MCLEAN, 

VA 

72102 

OICY 

ATTN 

J 

CCCKAYNF 

OICY 

ATTN 

B 

CHAMgfc  s  III 

OICY 

ATTN 

M 

KNASEL 

OICY 

ATTN 

W 

layscn 

OICY 

ATTN 

R 

SIFVERS 
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0F° AR  TMFN T  OF  OEFFNSF  CONTRACTORS 


SOlJTHWCST  RESEAPCF  INSTITUTE 
P  0  DRAWER  78510 
SAN  ANTONIO,  TX  78284 
OICY  ATTN  A  WFNZEL 
OICY  ATTN  W  Q  AKER 


SRI  INTERNATIONAL 
333  P  AVENSWOOf)  AVF'IUE 
MENLO  PARK,  CA  040?5 

OICY  ATTN  C,  ARRAH«WRCK 
OICY  ATTN  L  I  ?  P  A°  Y 
OICY  ATTN  J  CCLTCN 


SYSTEMS,  SCIENCE  S  SCETWARF  INC 
P  0  BOX  8243 
ALBUQUERQUE  *'M  °71«i8 
OICY  ATTN  C  NFFOHAV 


SYSTEMS,  SCIENCF  £  SrF  TWARE ,  INC. 
P  0  R^x  1620 
LA  JOLLA,  CA  07038 
OICY  ATTN  J  PARTf-cl 
OICY  ATTN  T  P  TNEY 
OICY  ATTN  I)  GRIME 
OICY  AT T N  LI  P ° A R Y 
OICY  ATTN  C  HASTING 
OICY  ATTN  K  PVATT 
OICY  ATTN  C  OTS''LXcS 
OICY  ATTN  T  CH^7> 


SYSTEMS,  SCIFNCU  f.  SOFTWARE  ,  INC. 
llQ00  SUNRISE  V A L L c Y  C^IVE 
RES  TOM,  VA  2  Zoc>  I 
OICY  ATTN  J 
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nco/iqTV'P^T  0C  DEFENSE  CONTRACTORS 


TEL  EDYME  BPOWN  P.NG!  NEF°  I\~. 
C'lMMTNGS  PESEA^CH  3APK 
HUNTSVILLE,  AL 

OICY  ATTN  J  R AVENSCp AFT 
OICY  ATTN  J  MfS.-.'Air 


TERRA  TEK,  INC. 


420  WAX AR A  WAY 
SALT  LAKE  CITY, 

LT  6 41 C  9 

OICY 

ATTM  A 

A  BO!  1- SAY  EC 

OICY 

ATTN  LIBRARY 

OICY 

ATTN  A 

JCN'E  S 

OICY 

AT^  S 

G0EEN 

TFT R A  TECH,  INC. 

430  M  POSEMEAO  Pl.VO 
° AS  AOFNA ,  CA  911^7 
OICY  A7TN  L  HV'ANC- 


TRW  OEFENSF  C  S°  AC.F  SYS  G°D'JP 
CINE  SPACE  PARK 
PEDONOn  BEACH,  CA  PC27  9 
OICY  ATTN  N  L  TPNER 

OICY  ATTN  TECHNICAL  IN^np^ATICM  CENTER 
OICY  ATTN  T  vaZZCLA 


TRW  DEFENSE  f.  SPACE  SYS  Gor,t'P 

p  o  «nx  mo 

SAN  BEPNARDIMC,  CA  <5240? 

OICY  ATTN  G  HVLCHFR 
OICY  ATTN  P  DAI 
OICY  ATTN  E  '-CNG 


C F  °A° r  ’■>r»,T  CF  OFF  c  NS  F  CCNT^ACTORS 


UNIVFPSAL  ANALYTICS,  IVC. 
7740  W  MANCHESTER  RlVC 
'’LAYA  OFl  REY,  CA  9C?cI 
0  ICY  ATTN  F  F  I E l.  0 


WEIOMNCEP  ASSOC..,  CCNSUITINC  CNGTMFFRS 
110  E  50 TH  STREET 
NFtf  YORK,  MY  100?’ 

01CY  ATTN  1  SANOLEP 

o ic y  attn  y  p  a r n n 


WF I  01  I  MGER  ASSOC.,  CC'SLLTIKC-  ENGINFFPS 
0000  SAND  HILL  °CAD 
MENLO  PARK,  CA  04075 
0 1C  Y  ATTN  J  ISENPEPO 


CJUjzJ 

'Tl&JidL  ']&ae<*suLk' 

$Cy<U*L<>> 

idcy  CUbC*.  ^ 
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0  C  P  A  R  T  M  F  NT  OF  OFFFMSF 

TO  TMF  CcCPFTAoy  r  r  OEFFNSC 
(  AT  CM  !C  ENERGY) 

WA  SHI  NOT*?*’  ,  DC 

OICY  ATTN  EXrCMTIVF  ASSISTANT 


OT SECTOR 

DEFENSE  com VI ink  ATTCNS  ACrMC v 
WASHINGTON.  DC  2’'',''c 

(A  r,o  c  ''V  ">  I  :  iU>  CrOE  ’40  FCR) 
OICY  A  T  TN  CnOF  S?n  C  ljdp 


DTC  EC  TOR 


OFFENSE 

INTEI  L  I  C. FMC  E 
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